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ABSTRACT 

We use the source counts measured with the Multiband Imaging Photometer for Spitzer (MIPS) 
at 24, 70, and 160 ^m to determine the 5-ct confusion Hmits due to extragalactic sources: 56/iJy, 3.2 
and 40 mJy at 24, 70 and 160 ^m, respectively. We also make predictions for confusion limits for a 
number of proposed far infrared missions of larger aperture (3.5 to 10m diameter). 
Subject headings: infrared: galaxies - galaxies: evolution - galaxies: statistics 



1. INTRODUCTION 

In addition to detector/photon noise, cosmological sur- 
veys in the far-infrared (FIR) spectral range are limited 
in depth by: 1.) structure in the infrared cirrus emission; 
and 2.) confusion by extragalactic sources. The first of 
these limitations can be avoided for some programs by 
observing in particularly low-background regions on the 
sky The second limitation arises because the high density 
of faint (resolved or unresolved) distant galaxies creates 
signal fluctuations in the telescope beam (Condon 1974; 
Franceschini et al. 1989; Helou & Beichman 1990; Rieke 
et al. 1995; Dole et al. 2003; Takeuchi & Ishii 2004, 
for instance). Because distant galaxies are distributed 
roughly isotropically and with a high density compared 
to the beam size, this noise is unavoidable. 

Extragalactic confusion noise can be robustly esti- 
mated by measurements of source counts combined with 
modeling to extend the counts to faint levels. We use 
new determinations of number counts in the three Multi- 
band Imaging Photometer for Spitzer (MIPS, Rickc ct 
al., 2004) band s, 24, 70, and 160 ^m (Dole ct al. 2004; 
iPapovich et al.ll2004t). and a m odel fitting all those ob- 
servables l|Lagache et al.ll2004D to determine more accu- 
rate limits for extragalactic confusion than have been 
available previously. Extragalactic confusion noise does 
not strictly follow Gaussian statistics. Therefore, we dis- 
cuss confusion limits in four different ways that are ap- 
propriate to various measurement situations: the pho- 
tomctr ic criterion and s ource density criterion (hereafter 
SDC) IjDole et al.l2003(l . and the levels deduced from the 
source densities of one source per 20 and 40 independent 
beams. We parameterize the noise as a "5 — cr" limit 
calculated as if it were Gaussian, because it is difficult to 
derive any other simple metric. All the definitions and 
values relative to MIPS Spitzer beams are summarized 
in Table 1 of Dole ct al. (200|). 
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We summarize the confusion limits for Spitzer in its 
three far-infrared bands in Tabled The situation is dif- 
ferent at 24 and 70 /im than at 160 /im. In the two 
first bands, where the background is resolved to a signif- 
icant extent, the confusion mainly results from the high 
density of resolved sources and their interference with 
extraction of fainter ones: the SDC is the appropriate 
measure (and the classical photometric criterion under- 
estimates the confusion level). In the third band, where 
the background is not well resolved, the confusion results 
from a population fainter than the sensitivity limit. In 
the latter case, confusion (and GIB fluctuation) proper- 
ties are directly linked to galaxy populations not directly 
detectable but which modulate the background level: the 
photometric criterion is appropriate. 

2. CONFUSION IN THE MID- AND FAR- INFRARED 

2.1. Confusion of Extragalactic Sources at 24 fxm 

The available measurements extend well into the ex- 
tragalactic confusion regime at 24/im, and the detector 
performance is also well understood even for long integra- 
tions. Therefore, we use this band to develop the general 
principles applicable to determining the confusion limits 
in Spitzer mid- and far-infrared imaging data. 

2.1.1. Confusion Limit Calculation 



Our confusion estimates a r e bas ed 
ogy described by iDole et all (120031^ . 
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number counts determined bv lPapovich et all (|200J) , ex- 
tr apolated to fainter flu x limits according to the model 
of iLagache et al.l Ij2004l) . Because these counts indicate 
that the background will be largely resolved into individ- 
ual sources, the appropriate measure of the confusion is 
the SDC. We obtain 56 /iJy for the 5-a confusion level, 
corresponding to 12 beams per source. It appears that 
this confusion level is in perfect agreeme nt with the 5a 
pre-launch predictions of lXu et al.l 1)200 Ij) . even if it was 
derived differently. If it were limited by photon noise 
only, the instrument would reach a detection limit of 
56 jiJ y 5-tT in 1900 seconds of integration (Rickc et aQ 
Hn^), so the model predicts that the gain in signal to 
noise ratio will have leveled out significantly for integra- 
tions of this length. 
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TABLE 1 
Spitzer MIPS Confusion Levels" 



time [s] 



Fig. 1. — Evolution of (Jtot (resulting contribution from the con- 
fusion noise and instrument noise, and derived from the Gaussian 
fit in the brightness map pixel histogram) as a function of integra- 
tion time, with a fit (dash) of the form: a?. = a^ . -I- cr-^ , „ = 

' V / tot znst conjHr 

At^^+C^. Dot-dash: constant term C Dot: \/A/t term. Top 
panel (a): 70 fim. Middle panel (b): 160 ^m. Notice the different 
scales in time (seconds) and atot (in brightness /i Jy/arcsec^ ) . 



There is excellent agreement between t he observed 80% 
completeness level and source density of IPapovich et all 
(|2b04i:) . and our SDC confusion level. However, it should 
be possible in principle to integrate below the 56 fiJy 
level, on a selected field of very low source density. In 
the "GOODS Test Field" in ELAIS Nl (described in Pa- 
povich et al. 2004), we estimate the area suitable for a 
deeper integration to be about 5% of the field area. 

2.L2. Noise Analysis 

We desired a test of these predictions that was as much 
as possible independent of assumptions about the in- 
frared galaxy population. For this purpose, we have char- 
acterized the noise in the 24 //m data from the ELAIS 
Nl field, the deepest observation obtained to date at this 
wavelength. We selected a very cleanly reduced region 
in the field, about 2x4 arcmin in size. We prepared 
two versions of the image in this region, both reduced 
identically, but one with an integration of 630 seconds 
and the other with an integration of 3800 seconds. We 
determined the pixel signal histogram in two ways. 1.) 
On a small region that also appeared to be free of de- 
tected sources, we verified that the standard deviation 
as measured in these histograms scaled inversely with 
the square root of the integration time. 2.) On the en- 
tire 2x4 arcmin region, we fitted it with a Gaussian 
of width fixed to the expectation for detector/photon 
noise. We required this Gaussian to fit only the nega- 
tive side of the histogram, on the assumption that there 
were no negative sources. We took the departure of the 
measured histogram from this fit toward positive fluc- 
tuations to be the infiuence of (at least) sources in the 
field. We measured the extension of the distribution to- 
ward positive values at half maximum. We found that 
the positive-going width of the distribution was larger 
than the pure detector/photon noise expectation by a 
factor of 1.7, in qualitative agreement with the effects of 
confusion. These excess fluctuations likely result from a 
combined effect of extragalactic sources, faint cirrus and 



24 /^m 



70 /^m 
[mjy] 



160 fj.m 

[mJy] 



SDC*" 


56 


3.2 


40 


20b<= 


71 


3.5 


45 


40b'* 


141 


6.3 


63 


Phot<= 


8 


0.7 


45 



Note. — (a) with Lajjachc ct al. 1 2004} model, (b) Using 
the Source Density criterion I Dole ct al. 2003') . (c) Using the flux 
corresponding to one source per 20 beams, (d) Using the flux 
corresponding to one source per 40 beams, (e) Using the standard 
photometric criterion and g = 4, for illustration. 



zodiacal light gradient. It is not clear at this stage which 
component dominates the fluctuations. 

2.1.3. Monte Carlo Simulation 

To empirically quantify the effect of confusion, we car- 
ried out a Monte Carlo simulation of source extraction 
under the conditions appropriate for the Spitzer deep 
24/im exposures. The approach is described in detail by 
Rieke et al. (1995). We built up a test field by distribut- 
ing confusing sources randomly according to a power law 
distribution matching the faint Spitzer number counts. 
Each source was entered as an Airy pattern. A test 
source of known amplitude was added to the center of 
the array, along with Gaussian noise. The sources were 
then identified using a modified GLEAN algorithm and 
finally the signal to noise was measured in a master ar- 
ray built up from the results of the GLEAN process, and 
in extraction apertures of various sizes. An important 
aspect of this simulation is that it combines the effects 
of neighboring bright sources and of the underlying, un- 
resolved distribution of faint ones, in a consistent man- 
ner. It should give a good measure of the confusion noise 
independent of the division between source density and 
photometric criteria. 

In the simulation, we excluded all objects brighter than 
400/iJy to avoid undo noise from bright-source artifacts. 
The first set of runs tested the extraction of a 56/iJy 
source in an 0.8 A/D beam, the beam size previously 
indicated to provide optimum performance in a heavily 
confusion-limited situation (Rieke et al. 1995 - this result 
was confirmed by the new calculations). We made 1200 
runs for an integration time long enough to drive detec- 
tor/photon noise down to 12.5/iJy, 5-ct. They yielded a 
net 5-cr limit of 60/iJy: removing the detector/photon 
noise leaves 59/xJy of confusion noise. That is, this ap- 
proach agrees well with the SDC-determined limit of 
56/iJy. 

We also simulated the results to be expected from 
shorter integration times. For example, if the 5-cr detec- 
tor/photon noise limit was set to 65yLtJy, then the indi- 
cated 5-cr level of confusion noise was 76/xJy, significantly 
poorer than from the simulation of very long integrations. 
This effect probably results from the increased uncer- 
tainty in source centroiding and the resulting lower ac- 
curacy in extracting accurate source measurements from 
a confused field. To test this hypothesis further, we sim- 
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ulated extraction of a 36/iJy source in the high signal-to- 
noise integration case, and found that the indicated 5-cr 
confusion hmit rose to 64/iJy, confirming the effect. 

2.2. Confusion by Extragalactic Sources at 70 jim 

At 70 /xm, we again use the number counts IjDole et al.l 
120041) as the basic input to determining th e conf usion 
level. The updated model of iLagache et alJ ^^2004^ was 
used to extrapolate the counts and to derive updated con- 
fusion limits. The use of a model is critical in this case 
because the contribution of unresolved sources is not neg- 
ligible. We derive a confusion level at 70 /im of 3.2 mJy 
using the SDC (Table ^. The differential source counts 
are almost flat (when divided by the Euclidean compo- 
nent), and the contribution from unresolved sources is 
much smaller than that of the resolved sources. These 
results demonstrate that the SDC estimate is the appro- 
priate one, that is, the confusion is dominated by faint 
resolved sources rather than the unresolved background 
due to even fainter objects. Further details are given in 
Tabled From the instrument radiometric model, we es- 
timate that about 1800 seconds of integration would be 
required to reach this limit. 

Again, we sought to check these results by a pure fluc- 
tuation analysis on the data without referring to galaxy 
population models. We used the data described by 
iDole et all ()2004(1 for the Chandra Deep Field South. We 
determined the evolution of Utot, the standard deviation 
of a Gaussian fltted to the surface brightness distribu- 
tion as a check of the results from extrapolating number 
counts downward. Data were combined into 6 mosaics 
corresponding to 100s to 600s integration time per sky 
pixel with 100s steps. Figure^ shows the evolution of 
CtotTO with time. We do not observe substantial flatten- 
ing in the crtono time evolution. We conclude that MIPS 
70 /im surveys do not reach yet the confusion limit after 
600s of integration. An estimate of the confusion level 
is given by fitting the time evolution of Utotro- We find 
that the detector/photon noise will be roughly equal to 
the confusion noise at > 800 seconds of integration, with 
large uncertainties because the fluctuation curve is still 
dropping almost as the square root of the integration 
time at the longest integration available. As at 24/im, 
this result is in satisfactory agreement with the integra- 
tion time predicted by the SDC modeling. 

2.3. Confusion by Extragalactic Sources at 160 fxm 

The data used at 160 ^m are also described by 
IDole et all l)2004j) . The ILagache et al . (2004) model pre- 
dicts a confusion level of 40 mJy (Table ^ . From the 
instrument radiometric model, we estimate that about 
70 seconds of integration would be required to reduce 
the instrument and photon noise to the level of the con- 
fusion noise. 

A similar fluctuation analysis as at 70 /im was con- 
ducted at 160 /im, where 6 mosaics corresponding to in- 
tegration times of 20 to 120s (with 20s steps) were stud- 
ied. Analyzing the fluctuations is more difflcult in this 
case since bright sources in the Euclidean regime con- 
taminate the statistics, and since the map S/N is not 
uniform. Nevertheless, we estimate from Figure Qd that 
the confusion and detector/photon noise should be equal 
at about 95 seconds of integration, in good agreement 
with the result from the SDC analysis. 




4 6 

Diameter [m] 

Fig. 2. — Confusion lev el vs Telescope Diameter, predicted by the 
So urce Density C r iterion IDol e ct al. 2003) with the updated model 
ofLa gache et alj 120041) . at 24 (plus), 70 (star) and 160 ^m (di- 
amond). Diameters refer to Spitzer, Herschel and SPICA, JWST 
and SAFIR. Dash: inverse square diameter law shown for illustra- 
tion. 



2.4. Confusion by Galactic Cirrus 

Another sensitivity limitation arises due to the struc- 
ture of the IR cirrus. To estimate how this cirrus emis- 
sion may affect the source detectability, we compared 
the 80% completeness limit in sky regions characterized 
by different cir rus background l evels, using simulations 
as described in iPapovich et al.l l|2004|) . We used a ded- 
icated engineering observation in Draco of a bright cir- 
rus, of HI column density nni varying between 4 and 
14 X lO^^cTO^^. At 24/im, we find a relatively weak effect, 
and derive a completeness degradation of 15% (~ 50/iJy 
increase from 340/iJy) between the dark and the bright 
parts of the cirrus field. The effects of cirrus are more 
conspicuous at 70/im. We reach a 80% completeness in 
Draco of '^17mJy and ~27mJy. In a low-cirrus field (e.g., 
Marano) and for a similar integration time (100s), this 
level drops at ^12mJy. We compared the estimates in 
Draco with those provided by the Performance Estima- 
tion Tool of the Spitzer Science Center and found that 
the measured value variations as a function of the cir- 
rus strength are in general agreement (within 30%) with 
those estimated by the tool from low to medium back- 
ground. This comparison will be refined as we continue 
to acquire far-infrared data. 

3. IMPLICATIONS FOR FUTURE OBSERVATORIES 

A number of cryogenically-cooled space telescopes have 
been proposed for the MIR, the FIR and the submil- 
limeter spectral ranges. Table |21 summarizes the main 
charac teristics of so me of the se observatorie s. Her- 
schel (IPilbratHl200l . JWST iCardned l200^. SPICA 
l)Matsumotoll2003() and SAFIR IjYorke et alJ l2002). have 
at least one photometric channel in common with MIPS. 
As examples, we focus on Herschel-PACS at 75 and 
170 /im, on JWST-MIRI at 24 /im, and on SPICA and 
SAFIR at 24, 70 and 160 /im, assuming in each case that 
the MIPS filters wifl be used. 

For each of these observatories, we compute predic- 
tio ns for the confusion le vel for unbiased surveys using 
the ILagache et al.l (2004) model of source counts. We 
assume a Gaussian beam profile for these future obser- 
vatories, with a FWHM of 1.22X/D, X being the wave- 
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TABLE 2 

Telescopes and Predicted Confusion Levels 



Hcrschcl" 
SPICA 



JWST'' 



SAFIR 



Diameter (m) 


3.5 


6.0 


10.0 


24 ^lm SDC ^Jy] 


2 


0.18 


<o.oi-f 


70 ^lm SDC mjy] 


0.16 


- 


0.004 


160 iim SDC* [mJy] 


10 


— 


0.6 



Note. — (a) With PACS. (b) With MIRI. (c) cf notes a and b 
1 Table HI (d) outside the range of the current model flux grid. 



In Table [SJ we use the confusion level given by the 
SDC, and compute the fraction of the CIB potentially 
resolved into sources. In the MIR, a significant step will 
be made with the 4ni-class space telescope: as an exam- 
ple, SPICA would potentially resolve 98% of the CIB at 
24 ^m. AU (> 99%) of the CIB would be resolved with 
JWST or SAFIR (although doing so with JWST would 
require extremely long integrations). In the FIR, Her- 
schel would resolve a significant fraction of the CIB at 
70 and 160 fim (resp. 93 and 58%, again with extremely 
long integrations). SAFIR will ultimately nearly resolve 
all of it (> 94%). 



TABLE 3 
Potential Resolution of the Cosmic Infrared Background 



24 ^m° 



70 /im" 



160 /xm" 



Spitzer 


74% 


59% 


18% 


HerschelVSPICA 


98% 


93% 


58% 


JWST*" 


99% 


- 


- 


SAFIR 


100% 


99% 


94% 



Note. — (a) Using the CIB value from Lagache et al. |2004) 
and using the limiting flux using the SDC limit and assuming 
confusion-limited surveys, (b) This hypothesis might not be valid 
for Herschel and JWST. 



length and D the diameter of the primary telescope mir- 
ror, given in Tabled The underlying assumption for the 
deepest surveys to be made by these planned facilities, 
is that they will be confusion-limited. This means we 
did not take into account other sources of noise, for in- 
stance photon noise due to insufficient integration times, 
or thermal background due to the warm telescope - by 
design, Herschel and JWST might be in the latter case. 
Normally background limited photon noise observations 
would give a sensitivity limit scaling as aperture squared 
for a diffraction limited system. Figure|21shows that con- 
fusion noise at 24 and 70 /im drops much faster than as 
aperture squared (dash) because source counts are shal- 
lower below fluxes where most of the CIB has been re- 
solved into sources. That is why the next generation of 
large far infrared telescopes will be much less confusion 
limited than Spitzer. 



4. CONCLUSIONS 

Using MIPS Spitzer data at 24, 70 and 160 um 
the source dens i ty me asured by iPapovich et alJ 1)20041) 
and iDole et alJ ll2004D together with the modeling of 
iLagache et al.l l^2004^ has allowed us to derive the con- 
fusion limits for Spitzer in the mid to far infrared. We 
tested the model results with a Monte Carlo simulation 
at 24/zm and with a fluctuations analysis at all 3 wave- 
lengths. The agreement is uniformly very good. 

At 24 and 70 /Ltm, confusion is mostly due to the high 
density of resolved sources, and at 160 /im, confusion is 
mainly due to faint unresolved sources. Studying the FIR 
fluctuations at this wavelength is thus a tool to constrain 
the nature of the faint galaxies, beyond the confusion 
limit. 

We also derive confusion limits for future space IR 
observatories. We show that future large-aperture mis- 
sions will gain in confusion-limited sensitivity substan- 
tially faster than as aperture squared for wavelengths < 
100/im, allowing them to reach very deep detection lim- 
its. For example, the CIB should be fully resolved into 
sources in the MIR and FIR with SAFIR observations. 
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